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Abstract. BPS invariants are computed, capturing topological invariants of moduli spaces 
of semi-stable sheaves on rational surfaces. For a suitable stability condition, it is proposed 
that the generating function of BPS invariants of a Hirzebruch surface takes the form of a 
product formula. BPS invariants for other stability conditions and other rational surfaces are 
obtained using Harder-Narasimhan filtrations and the blow-up formula. Explicit expressions 
are given for rank < 3 sheaves on and the projective plane P'^. The applied techniques 
can be applied iteratively to compute invariants for higher rank. 



1. Introduction 

Topological invariants of moduli spaces of semi-stable sheaves on complex surfaces are a 
rich subject with links to many topics in physics and mathematics. Closely related topics in 
physics are gauge theory, instantons, electric-magnetic duality [26] and also (multi-center) 
black holes [H HSl ISO]- Instantons saturate the bound on their minimal action, the so-called 
Bogomolnyi-Prasad-Sommerfeld (BPS) bound. The prime interest of this article are topo- 
logical invariants of moduli spaces of instantons, in particular their Poincare polynomials, 
which are commonly referred to as "BPS invariants". These invariants correspond also to 
(refined) supersymmetric indices enumerating supersymmetric or BPS states. 

Instantons on complex surfaces correspond algebraically to semi-stable vector bundles and 
coherent sheaves [TT| |5]. Generating functions of BPS invariants of sheaves on surfaces are 
computed for rank 1 by Gottsche [7j and rank 2 by Yoshioka fIE\ 12^) . These generating 
functions lead to intriguing connections with (mock) modular forms PB| El El 12], which are 
a manifestation of electric-magnetic duality of the gauge theory f2U\. Refs. |2H 123] compute 
BPS invariants for rank 3 sheaves with Chern classes such that stability coincides with 
semi-stability. 



2000 Mathematics Subject Classification. 14J60, 14D21, 14N35. 
Key words and phrases, sheaves, moduli spaces. 

1 



2 JAN MANSCHOT 

The present article computes the BPS invariants of semi-stable sheaves with rank 3 on 
Hirzebruch surfaces and on the projective plane P^, and explains how to generalize the 
computations to higher rank. The developed techniques can be applied straightforwardly to 
compute BPS invariants of the other rational surfaces. Although the extension from stable to 
semi-stable might seem a minor one, it requires to deal with various subtle but fundamental 
aspects of the moduli spaces of semi-stable sheaves, which could be neglected in Ref. |21] . 
Having resolved how to deal with these aspects for r = 3, the computations can in principle 
be extended to any rank. 

This introduction continues with summarizing the contents of the paper, after recalling 
the computations in Ref. |21] which were inspired by [28| [29| [H [9j. A crucial fact for the 
computations is that the blow-up : — )■ is equal to the Hirzebruch surface — )■ (7 
with i = 1. The fibre / and base C of are both isomorphic to P^. Since no stable bundles 
on P^ exist, the BPS invariants of with polarization J chosen sufficiently close to / (a 
so-called "suitable" polarization, see Definition 15. ip vanish for sheaves with first Chern class 
Ci and rank r such that Ci ■ / 7^ mod r. 

Wall-crossing then allowed to compute the BPS invariants for other choices of J. The 
BPS invariants of P^ were obtained from those of P^ by application of the blow-up formula 
[30| [9l [18], which is a simple relation between the generating functions of the invariants for 
P^ and P^. However, its original form is only valid for gcd(ci ■ (f)*H,r) = 1 and J = (l)*H, 
with H the hyperplane class of P^. 

The present paper describes how to deal with the cases when Ci and r do not satisfy the 
constraints for vanishing of the BPS invariant or the blow-up formula. The formal theory 
of invariants of moduli spaces (or stacks) of semi-stable sheaves is developed by Joyce [12j. 
We will in particular use the notion of Poincare functions for moduli stacks, which are a 
generalization of Poincare polynomials of manifolds. The Poincare function of a moduli 
stack is (conjecturally) related to the BPS invariant by fl3.5p . The BPS invariant is most 
natural from physics and leads to generating functions with modular properties. 

The two novel ingredients of this paper are: 

(1) Eq. (14.21) which provides for any rank r > 1 the generating function of Poincare 
functions of the moduli stack of sheaves on P^ whose restriction to the fibre / is 
semi-stable. 
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(2) Extended Harder-Narasimhan filtrations C -Fi C -F2 C ■ • • C -F^ = -F, whose 
definition (Def. I5.3p differs from the usual definition (15.21) of HN filtrations by allowing 
quotients Ei = Fi/Fi^i with equal (Gieseker) stability n) ^ pj{EiJ^i,n). These 

filtrations in combination with the associated invariants fl5.8p are particularly useful 
to compute generating functions of BPS invariants starting from Conjecture 14.21 and 
their changes accross walls of marginal stability. 

To obtain the BPS invariants for a suitable polarization, one subtracts from Eq. (14.21) 
generating functions corresponding to extended HN-filtrations given by (15. Sp . analogous to 
the seminal papers about vector bundles on curves [T0| [1]. Naturally, these techniques are 
also applicable to compute invariants of semi-stable invariants for other mathematical objects 
like vector bundles on curves and quivers. Also a solution to this recursive procedure is given 
analogous to Ref. [52]. Then repeated application of the formula for filtrations (which is 
equivalent with the wall-crossing formulas |16 [ [T3 ] ) gives the BPS invariants for other choices 
of the polarization. 

Finally, the blow-up formula provides the invariants on P^. The earlier mentioned condi- 
tion gcd(ci ■ (f)*H, r) = 1 is a consequence of the fact that the blow-up formula is applicable 
for the Poincare functions T^{r, w; J) with respect to /i-stability instead of the more refined 
Gieseker stability. However with the invariant for filtrations (15. 8p . it is straightforward to 
transform the BPS invariants Q{r,w;J) to X^(r,w; J) for /i-stability. The rational fac- 
tors in Eq. (15.81) appear naturally in the relation between the generating functions of these 
invariants. 

The paper illustrates in detail the above steps for sheaves with rank 2 and 3, and shows 
their agreement with various consistency conditions, e.g. the blow-up formula, integrality 
and w ^ symmetry of the Poincare polynomial. 

The outline of the paper is as follows. Section [2] reviews some necessary properties of sheaves 
on surfaces including stability conditions. Section [3] defines the invariants and generating 
functions. Section H] presents the formula (14. 2 p for the contribution to the BPS invariants 
from sheaves whose restriction to the fibre is semi-stable. Then we continue with the com- 
putation of the invariants of for any choice of polarization in Section [5l To this end an 
invariant is associated to a filtration with prescribed topological classes of the quotients. 
Finally Section [6] computes the generating function for sheaves on P^ with (r, Ci) = (3,0). 
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2. Sheaves on surfaces 

We consider sheaves on a smooth projective surface S. The Chern character of the sheaf 
F is given by ch(F) = r{F) + Ci(-F) + ici(-F)^ — C2{F) in terms of the rank r{F) and its 
Chern classes Ci{F) and C2{F) . The vector T{F) parametrizes in the following the topological 
classes of the sheaf T{F) := ( r{F), chi(F), ch2(-F) ). Other frequently occuring quantities are 
the determinant A(F) = ^{c2{F) - r^ci{Ff), and fi{F) = ci{F)/r{F) G H\S,Q)- 

Given a filtration C Fi C ■ ■ ■ C = F, let = Fi/Fi_i and Ti = T{Ei). The 
discriminant of F is given in terms of the subobjects and quotients by: 

P.l) A( r(F) ) = g ^A(S) - 5^ t '"'^'('g-'' (P(F.) - A(fi-O)^ . 

We are interested in the moduli space (or moduli stack) of semi-stable sheaves with respect 
to Gieseker stability, but also the coarser ^u-stability appears in order to apply the blow-up 
formula. To define these two stability conditions, let C{S) C H'^{S,M.) be the ample cone of 
5", and the (reduced) Hilbert polynomial pj{F,n) = x(F ® J"-)/r{F). For a surface S, we 
have [5]: 
(2.2) 

X .9 9/ fci(F)-J Ks-J\ 1 f cAF)^ - Ks ■ ci(F) , , 
p,(F, n) = J^nV2+ ^-^^ ^ j ( ^ ^ 2 " ''^^^) ^^^^^^^ 

Note that this function can be obtained from the physical central charge as in [H In 
the large volume limit, the stability condition asymptotes to the lexicographic ordering of 
polynomials based on their coefficients. This ordering is denoted by -<. Then, 

Definition 2.1. A torsion free sheaf F is Gieseker stable (respectively semi-stable) if for 
every subsheaf F' C F, pj{F',n) -<pj{F,n) ( respectively pj{F' ,n) :<pj{F,n) ). 



and 
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Definition 2.2. Given a choice J G C{S), a torsion free sheaf F is called ^-stable if for 
every suhsheaf F' C F, fi{F') ■ J < fi{F) ■ J, and fi-semi-stable if for every sub sheaf F' , 
jj{F') ■ J < 12(F) ■ J. 

Thus /^-stability is a coarser stability condition then Gieseker stabihty, although the walls 
for both stability conditions are the same. We will find in Sections |5] and |6] that invariants 
based on Gieseker stability exhibit better integrality and polynomial properties then the ones 
based on /i-stability. On the other hand, operations like restriction to a curve and blowing- 
up a point of S are most natural for /i-semi-stable sheaves. A wall of marginal stability 
W{F', F) C H^{S, M) is the codimension 1 subspace of C(5), such that {^{F')- ^{F))- J = 0, 
but (/i(F') - • J ^ away from W{F', F). 

The moduli space Aij{T) of Gieseker semi-stable sheaves on 5* (with respect to the ample 
class J) whose rank and Chern classes are determined by F has complex dimension: 

(2.3) dime Mj{T) = 2r^A - r\{Os) + 1. 

Twisting a sheaf E hy a line bundle C is an isomorphism of moduli spaces. The Chern 
classes of the twisted sheaf E' = E ® C are: 



r{E') = r{E), ci(E') = ci(E) + r(E)ci(£) 
C2{E') = C2{E) + {r{E) - l)ci(£)ci(E) + ci(£) 



2r(E)(r(E)-l) 



2 

The discriminant remains invariant: A{E') = A{E) . This shows that it suffices to compute 
the generating functions for Ci{E) G H"^ {S , Z / rZ) . 

Determination of the generating function for r > 2 is complicated in general. To make 
progress, we specialize in the following to the set of ruled surfaces. A ruled surface is a 
surface "Eg^i together with a surjective morphism vr : — )• C to a curve C with genus g, 
such that the generic fibre is isomorphic to for every point of C. Let / be the fibre of 
TT, then 7/2(^3/, Z) = ZC © Z/, with intersection numbers = —i, = and C ■ f = 1. 
The canonical class is K-^^^ = —2C + {2g — 2 — i)f. The holomorphic Euler characteristic 
x{^T.g e) is 1 — f?. An ample divisor J G C(Sp^£) is parametrized by J^^ = m{C + if) + nf 
with m,n > 0. 

This article will further specialize to the Hirzebruch surfaces Sq.^ = S^, and in particular 
to Si. Besides being a ruled surface. Si is also the blow-up : — )■ of the projective 
plane P^. The exceptional divisor of (p is C, and the pullback of hyperplane class of 
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equals C + f. Due to the simplicity of P^, it is of intrinsic interest to determine the 
generating functions of its BPS invariants. 



3. BPS INVARIANTS AND GENERATING FUNCTIONS 



This section defines the generating functions of the BPS invariants and discusses some 
of its properties. Although this article computes the BPS invariants and shows that they 
satisfy many desired and expected properties, their precise definition is still not completely 
settled. Physically the BPS invariant is given by a weighted sum over the BPS Hilbert space 
T-L^r, J). We choose to divide by a factor w — w^^ which appears physically as an overall sum 
of oscillators. The BPS invariant arises by considering topologically twisted A/" = 4 Yang- 
Mills on the surface 5* f2U\. The path integral of this theory localizes on the BPS solutions, 
in particular instantons, due to the topologically twisted supersymmetry [26j. The path 
integral shows in particular that the (numerical) BPS invariant corresponds to the Euler 
number of the BPS moduli space. 

Alternatively one can consider the M = 2 supersymmetric theory in M^'^ obtained from 
the compactification of IIA theory on the non-compact Calabi-Yau 0{Ks) — ?■ S. From this 
last perspective the BPS invariant can be refined with an additional parameter w j6]: 



with J3 a generator of the SU{2) = Spin(3) group arising from rotations in M^'^, and 
J3 is a generator of the S't/(2)/j /^-symmetry group. BPS representations have the form 



(the half-hypermultiplet) present for every BPS state [6j. Since Q{r,w; J) is thus essentially 
an SU{2) character, this shows that Q(r,w; J) is a polynomial divided hj w — w~^; the 
polynomial has integer coefficients and is invariant under w ^ w~^. The positivity conjec- 
tures of Ref. [6] assert furthermore that the coefficients are positive. Many formulas take a 
simpler form in terms of the rational BPS invariants given by: 



(3.1) 



n(T,w;J) 




[(|, 0) © (0, |)] ®a; with cu = a vacuum representation of Spin(3) (BSU{2)r with spins 

j and j'. One factor of w — w~^ in the denominator will vanish due to the factor (|, 0) © (0, |) 



(3.2) 




m|r 



They appear in fact rather natural from the physical perspective [22| [H]. 
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The M = 2 picture shows that the refined BPS invariant provides more information than 
the Euler number of the moduh space. The i(7-expansion is expected to give the Poincare 
polynomial of the BPS moduli space. To make this more precise, we let M.j(T) be the 
suitably compactified moduli space of semi-stable sheaves on 5* with topological classes F 
and for polarization J G C(S'), i.e. the Gieseker-Maruyama compactification. If semi-stable 
is equivalent to stable the moduli space is smooth and one can take as mathematical definition 
of the BPS invariant 

-dime A4j(r) 

(3.3) n{T,w-J) := —p{Mj{T),w), 

w — 

where p{X,w) = Yli=o^^^^^ hiX) is the Poincare polynomial of X, with the Betti 

numbers &j(X) = dimi7*(X, Z). Eq. (12. 3p provides us with the degree oi p{M.j{T),w), and 
since AijiV) is compact, orientable and without boundary hi{X) = b2dimc{x)-i{X). 

If semi-stable is not equivalent to stable, A4j(r) contains singularities due to non-trivial 
automorphisms of the sheaves, and the mathematical definition of the BPS invariant is not 
clear. The proper way to deal with the singularities is to consider instead of the moduli 
space A^j(r) the moduli stack VJtjiT). Let X{r,w; J) be the virtual Poincare function of 
the stack DJlj{r,w) of semi-stable sheaves \12\ (see also |16[ I2l]). The virtual Poincare 
function X{r,w; J) is a rational function in w and a natural generalization of the Poincare 
polynomial of smooth projective varieties to stacks. 

The seminal papers [U EH |29] compute moduli space invariants by explicitly counting 
sheaves on the surface S defined over a finite field with s elements. The Poincare function 
X(r, s^; J) is upto an overall monomial computed by: 

1 



(3-4) E 



, ,#Aut(E)' 

where Mj{r,¥s) is the set of semi-stable sheaves with characteristic classes F. The Weil 
conjectures imply that the expansion coefficients in s are the Betti numbers of the moduli 
spaces. The parameters s is related to the w in this paper hj s = w"^. Eq. (13. 4p shows that 
poles of I{r,w]J) in w appear when the sheaves have non-trivial automorphism groups. 
If semi-stable is equivalent to stable X(F,w; J) = Q{r,w; J); the factor w — is due to 
the automorphisms which are multiplication by C*. The automorphism group of semi-stable 
and unstable bundles or sheaves is in general GL{n), whose number of elements over is 
(1 — s)(l — s^) . . . (1 — s") and thus lead to higher order poles. 
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If semi-stable is not equivalent to stable, the following assumes (and gives evidence) that 
X(r, w\ J) is related to the multi-cover invariants VL(T , w; J) by: 

1 ^ 

(3.5) Ar,w;j)= Yl j^_l[n{r,,w;j), 

ri + ---+r^=r ' i=l 

pj(rj,n)=p j(r,n) for i = l,...,e 

with inverse relation: 

(3.6) Qir,,w;J)= Yl ^-^l[liT,,w;J), 

ri + .-.+r^=r i=i 

Note that X(r,t/;"^; J) ^ — X(r,w; J) and that X(r,ty; J) might have higher order poles in 
w as discussed above. The invariants X(r,w; J) are not so natural from the perspective of 
modular generating functions of BPS invariants, since the generating functions of r2(r, w\ J) 
have weight —h2{S)/2. However, the product formula of Conjecture 14.11 is a generating 
function for X(r, w\ f) and invariants X(r, w\ J) with respect to /i-stability instead of Gieseker 
stability appear in the blow-up formula in Section O 

Since the invariants X(r,w; J) are not so natural for generating functions, the following 
will mainly work with the invariants fi(r, J) and VL{T, w; J). This requires to multiply by 
factorials due to Eq. (13. 5p . Interestingly, these naturally combine generating functions of 
Gieseker unstable and strictly semi-stable sheaves. 

We will find that Q{r,w;J) indeed has the form (13. 3p . with the expected integrality 
properties. This generalizes Refs. |2S] where it is shown that Cl{r,w;J) computes 

the intersection cohomology of the singular moduli space Aij{r) for r = 2. Intersection 
cohomology is a cohomology theory for manifolds with singularities which satisfies Poincare 
duality if the manifolds are complex and compact. It is therefore natural to expect that 
the BPS invariants for r > 3 also provide Betti numbers of intersection cohomology groups. 
This issue is left for further research. 

We continue now by defining the generating function hr^ci{z, t; S, J) of Cl{r, w; J): 

(3.7) K,cA^, T- S,J) = J2 ^(r, w; J) g'^^(r)-^. 

where q := e^'^*'^, with r G and w := e^'^*^ with z & C Since twisting by a line bundle 
(12. 4p is an isomorphism of moduli spaces, it suffices to compute hr^ci{z,T; S, J) for Ci G 
H2{S, Z/rZ). The expansion parameter t for C2 in Refs. jZl EH |29] is related to g by g = s'^t. 
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The generating function /ii ^(^z, t; S*) depends only on b2{S) for S a smooth projective 
surface with bi{S) = 63(5*) = [7J: 

(3-8) /^i.c,(^,r;S) = ^-^^—^-^-^^-^ ^ 

where the Dedekind eta function rjlr) and Jacobi theta function 9i{z, r) are defined by: 

00 

Vir) := q^^U^l-qn, 

n=l 

e^{z,T) := iq^w'^ -w-^)Y[{l-q''){l-wq''){l-w-^q''). 

n>l 

The dependence on J is omitted in Eq. (13. 8p . since all rank 1 torsion free sheaves are 
stable throughout C{S). Similarly, J is omitted in the following from hr^ci{z, t; P^, J), since 
62 (P^) = 1 and therefore the BPS invariants do not vary as function of J. For clarity of 
exposition, is omitted from the arguments of hr,ci{z,T; S^, J). 

4. RESTRICTION TO THE FIBRE OF HiRZEBRUCH SURFACES 

This subsection deals with the set Mj(r) of sheaves whose restriction to the (generic) fibre 
/ of vr : S£ — )■ C is semi-stable. Since / = all vector bundles on / are isomorphic to a sum 
of line bundles, and therefore the restriction of a sheaf to the fibre can only be semi-stable 
if its degree d{E\f) = Ci ■ f equals mod rl3 For these Chern classes the Q(T; w, J) for 
J = Jo,i or for a suitable polarization J^ i (see Definition 15. ip . 

The BPS invariants for J^^i can be computed from the Poincare function X{r,w;f) of 
the moduli stack 971/ (F) of sheaves whose restriction to / is semi-stable. We define the 
generating function Hr{z, r; /) of X(F, w; f) by: 

(4.1) i/,(z,r;/):=5^X(F,t.;/)g'^^(r)- 

The following conjecture gives Hr{z, r; /) for any r > 1: 
Conjecture 4.1. The function Hr{z,T;f) is given by: 

ii-irzivirf 



X(S) 
24 



(4.2) Hr{z,T; f) = { ei(2^,r)2ei(4^,r)2...0i((2r-2)^,r)2ei(2r^,r)' if Ci ■ / - mod r, r > 1 

0, if Ci ■ / 7^ mod r, r > 1. 



Recall that a vector bundle F of rank r and degree d on a curve C is stable (respectively semi-stable) if 
for every subbundle F' C. F (with rank r' and degree d') d! /r' < d/r (respectively d' /r' < d/r). 
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Conjecture 14.11 is obviously in agreement with Eq. (13 .Sp for r = 1. Ref. [29j proved 
the conjecture for (r, Ci) = (2,/), which is now briefly recalled. Ref. [29j considers the 
ruled surface over a finite field F^, and utilizes the fact that any vector bundle in F 
can be obtained from 7r*7r*F, which is a vector bundle on supported on C, by successive 
elementary transformations. 

An elementary transformation is defined by 

Definition 4.2. Let D he an effective divisor on the surface X . If F and G are vector 
bundles on S and D respectively, then a vector bundle F' on S is obtained by an elementary 
transformation of F along G if there exists an exact sequence: 

(4.3) ^ F' ^ F ^ i^G ^ 0, 
where i denotes the embedding D G S . 

This shows that the contribution to h2^ci{z,T] J) from Mj(r) is the product of the the 
total set of vector bundles on G, multiplied by the number of elementary transformations. 
The total set of vector bundles with r = 2 on C is enumerated by: 

(4.4) f^Cc{2) 

1 — s 

with Cc{^) the zeta function of the Riemann surface G: 

(4.5) Cc{n)- ^^'-'^ 



[1 - s-")(i - s^-'^y 

Multiplication by the factor due to elementary transformations gives 

(4.6) E ^A^m = r=7^^(') n ^«(^' s'^-H'^)z.{s^ s^^n 

C2 EeMf{2,mf,C2) ^ ' a>l 

with Zs{S,t) the zeta function of the surface S: 

(4.7) Zs{S,t) 



(l-t)(l-St)''2(S)(l_s2i)- 

The parameter substitutions q = s^t and w'^ = s give then Eq. fl4.2p (upto an overal 
monomial in w and q). It is conceivable that Conjecture 14.11 for r > 2 can be proven in a 
similar manner. The following sections show that at least for r = 3, 4, it is consistent with 
various other results. Moreover, it continues to hold for the other Hirzebruch surfaces with 
£ > 0. 
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As an aside, we mention that the Poincare function of the stack of vector bundles on a 
curve Cg with genus g is given in terms of w by |T0| [1] : 

(4.8) H,{z,Cg).- w ^_ II (i_^2,)2 • 



The first term in the g-expansion of Eq. (14. 2 p starts with Eq. (14. 8p for g = 0. In fact, 
one could view Hr{z,T;f) as an extension of Hr{z;Cg) to a modular infinite product. It 
is conceivable that the generalization of Eq. (14. 2 p to ruled surfaces with g > is similarly 
related to Hr{z; Cg). For example, this gives for r = 1: 

(4-9) /.,,(.,r;S,,) = ^^^2,,.),(.)i.2.- 

5. BPS INVARIANTS OF HiRZEBRUCH SURFACES 

5.1. BPS invariants for a suitable polarization. This subsection computes for ci ■ / = 
mod r the BPS invariants of for a polarization J G C(S^) sufficiently close to J^ i = f. 
The BPS invariants are for this choice of J independent of i. "Sufficiently close" depends on 
the topological classes of the sheaf; it is known as a P-suitable polarization: 

Definition 5.1. A polarization J is called T-suitable if and only if: 

- J does not lie on a wall for T = (r, chi, ch2) and, 

- for any semi-stable subsheaf F' G F with T{F) = T, {fi{F') — /i(-F)) ■ / = or 
(/i(-F') — /u(F)) ■ / and {^{F') — /u(-F)) • J have the same sign. 

We will keep the dependence on the Chern classes implicit in the following and denote a 
suitable polarization by J^^i with e positive but sufficienty small. From the definition follows 
that if Je^i is an r(F)-suitable polarization, and F\f is unstable, then F is /i-unstable. Thus 
we need to subtract from M/(r), i.e. the set of sheaves with topological classes P whose 
restriction to the fibre / is semi-stable, the set of sheaves which are semi-stable for J = / 
but Gieseker unstable for i. We continue by explaining this for r = 2. Then the general 
formula is proposed for the invariant enumerating extended HN-filtrations in a suitable way, 
which is then illustrated for r = 3. 

A crucial tool to obtain the invariants enumerating semi-stable sheaves are Harder-Narasimhan 
filtrations [lOj, which can be defined for either Gieseker or yU-stability. To define these fil- 
trations, let (f denote either Gieseker, <f{F) = pj{F,n), or /^-stability, ip{F) = fi{F) ■ J. 
Then: 
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Definition 5.2. A Harder-Narasimhan filtration (HN- filtration) with respect to the stability 
condition if is a filtration G Fi G F2 G ■ ■ ■ G Fi = F of the sheaf F such that the quotients 
Ei = Fi/Fi-i are semi-stable with respect to (f and satisfy f{Ei) > ip{EiJ^i) for all i. 

Since /x-stability is coarser then Gieseker stability, the length Ig{F) of the HN-filtration 
with respect to Gieseker stability is in general larger than the length ifj,{F) of its HN-filtration 
with respect to //-stability. 

The Poincare function of the stack of HN-filtrations with respect to Gieseker stability and 
prescribed Fj = T{Ei) is [30j : 

e 

(5.1) w-^'<^'^''^^f'^-^'^^-^^Yll{r,w; J), 

1=1 

where rirj{fij — fii) ■ Ks is the Euler form for semi-stable sheaves on the projective surface 
S. One could define a similar function for the stack of filtrations with respect to yU-stability. 

For (r, Ci) = (2, /), the only HN-filtrations have length = 2. Denoting Ci{E2) = bC~af, 
and thus Ci{Ei) = —bC + (a + 1)/, one easily verifies that the HN-filtrations correspond to 
a > and 6 = 0. Since 6 = the dependence of Ks in Eq. (15.11) does not lead to a 
dependence on i. Using that Eq. (13.81) is also the generating function of X(F, w; J) for r = 1, 
Eq. (15. ip becomes: 

(5.2) 5^^-^(^'^+^)/^i,o(z,r)2 = --^^h,,o{z,rr, 

a>0 

where we assumed \w\ > 1. Subtracting this from Eq. (14. 2 p for r = 2 gives: 
(5.3) 



which is easily verified to enumerate invariants ^(F; J^^i) satisfying the expected properties 
mentioned below Eq. (13. 3p . 

For (r, Ci) = (2, 0), the HN-filtrations with = 2 split naturally in two subsets: the first 
set has length = 2 with respect to /i-stability, and the second set has = 1. Similarly to 
(15. 2p . the first set gives rise to: 



(5.4) -^—-^h,,o{z,Tr 
and the second set to: 



(5.5) h^,oiz, r)'-lYl ^((1' 0' 



n>0 
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where the second term subtracts from the first the Gieseker semi-stable sheaves which should 
not be subtracted from H2{z, r; /). Subtraction of Eqs. (15 ■4p and (I5.5p from H2{z, r; /) gives 
the generating function of X( (2, 0, C2), w; J), which corresponds by Eq. (13. 5p to: 

^^■^^ ^2,0(2;, t; Je 1) = —— ( /J ?! N + 1 ~ 7; 

0i[2z,T)'' ri[T)'^ \0i[4:Z,T) 1 — 2^ 

Again one can verify that the invariants satisfy the expected integrality properties. Ref. [29\ 

computes the Betti numbers of the intersection cohomology of the singular moduli spaces 

and arrives at the same generating function (15. 6p . 

The Betti numbers for the intersection cohomology of the moduli space of semi-stable 

vector bundles on Riemann surfaces were earlier computed in Ref. |15| . The above procedure 

gives these Betti numbers with much less effort. For example, one can easily verify that 

(5.7) H^iz, Cg) + (-^ - l) H,iz, Cg)\ 



l-w^ 2^ 

with Hr{z, Cg) as in Eq. (14. Sp . is equivalent with Proposition 5.9 in the Corrigendum to [15] . 

Since the invariants X(F,w; J) are not so compatible with modular generating functions 
for r > 2, it is useful to work as much as possible with the invariants VL{T , w; J). To this end 
an extension of the HN-filtration is necessary: 

Definition 5.3. j4n extended Harder-Narasimhan filtration (with respect to Gieseker stabil- 
ity) is a filtration G Fi G F2 G ■ ■ ■ G Fi = F whose quotients Ei = Fi/Fi_i are semi-stable 
and satisfy p J {Ei,n) >z pj{Ei+i,n). 

From Eq. (13. 5p follows that the natural invariant ^({Fj},^; J) associated to the stack 
97t({Fj}, ty) of extended HN-filtrations with prescribed Chern classes Fj = F(i?j) is: 

1 ^ 
(5.8) n{{V,}-w,J) := |Aut({Fj; j^i ^'^'^-^'^^^"^'"'^''"^ H ^(^^ J)- 

The number |Aut({Fj}; J)| is equal to Ha^^a'' where is the total number of quo- 
tients Ei with equal reduced Hilbert polynomial pj{Ea,n). Thus only for HN-filtrations 
|Aut({F,};J)| = l. 

If the sum over all extended HN-filtrations contains a group Ei with equal pj{Ei,n) but 
unequal Fj, the factor |Aut({r }-j)| divides out a number of permutations. To avoid this 
overcounting, one could introduce a further ordering on the vectors Fj, which should be 
obeyed by the set of filtrations to be summed over. Then one would divide by |Aut({Fj})| = 
J^pTip!, where rip is the number of equal vectors Fp appearing among the Fj, i = 
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This is the origin of the "Boltzmann statistics" in wall-crossing formulas [22j in the work of 
Joyce [12] . 

The functions /ir,ci(-2, t'-, Je,i) with Ci ■ / = mod r are given by the recursive formula 



rx(.S) 
24 



(5.9) K,c, {z, T- J,,i) = Hr{z, T-f)-J2 E ^(^rj; w, J,,i) g'-^(r) 

ch2 ri + ...+r^=(r,ci,ch2) 

Pj(ri,n)hPj(ri + l,n),t>l 

with A(r) given in terms of Fj by Eq. (12.11) . 

We continue by applying Eq. (15.81) to compute h^ci{z,T] Ji,e), with Ci = / and 0. One 
obtains: 

Proposition 5.4. 

(5.10) hJiz,T■J,,^) = ^^2^/l' ^2^... + 



0i(22,r)2^i(4^,r)2^i(6z,r) 1 - ^6 ^i(2z, r)3 ^i(4z, r) 

4 



(5.11) A3.(..r;4o ^ .... , 



^i(22,r)2^i(4z,r)2 0i(6;z,r) 1 - ^i(2z, r)3 ^(4z, r) 
1' 



(l-«;4)2 3; ^i(2z,r)3r7(r)3' 

Proof. We start by proving Eq. (I5.10p . Denote the length of an extended HN- filtration by 
i, its length with respect to /i-stability by and Gieseker stability iQ. We first consider 
the unstable filtrations with i = = 2, and parametrize Ci{E2) by bC — af. These are 
parametrized by a > and 6 = 0. There are four possibilities to be distinguished: whether 
r{Ei) = 1 or 2, and whether the quotient with rank 2 has Ci = or / mod 2. Adding up 
these contributions, one obtains: 

w'^ + up' + llP^ 

(5.12) - — hi,o{z, r) h2,o{z, r; J^.i) - t^t hi,Q{z, r) h2j{z, r; J^^i), 

1 — 1 — 

The filtrations with i = 3 consist of 3 subsets: one set with = 3, one with £^ = 2 but 
ic = 3, and one with iQ = 2. Parametrizing Ci(-E'j) = biC — aif, the first set is parametrized 
by tti — Cj+i > 0, Yl^=i (^i = ^ ciiid bi = 0. These are counted by: 

(5.13) Yl w-'^'^^''^^ h,o{z, rf = ^ h,o{z, rf. 

ki.k2>0 ^ ' 

fc2 — fc]^ — 1 mod 3 
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For the second and third sets, one needs to distinguish between equahty of the stabihty 
condition of E2 with Ei or E3. These two sets are enumerated by: 

(5.14) .^.'^h^„^,,rf. 

Note that the factor |Aut({r^} Ji )| naturally combines the contributions of filtrations with 
ifj, < i. Another observation is that the term — | in the second factor of /i2,o(-2, t; J^^i) (15. 6p 
cancels against (15 .Mp in the total sum. After subtraction of the terms (I5.12p -f l5.14p from 
Eq. (14. 2 p for r = 3, and writing the whole series in terms of modular functions, one obtains 

For (r, Ci) = (3, 0), one needs to subtract the following terms: 
- due to unstable filtrations with £ = £^ = 2: 

2 2w^ 



1 — 1 — w^^ 

due to unstable filtrations with £ = 2, £^ = 1 and £q, = 1 01 2: 

2 

due to unstable filtrations with £ = £^ = 3: 

hi,o{z,r) , 



- due to unstable filtrations with £ = 3, = 2 and = 2 or 3: 

~ o 1 12 ^1.0(^5 

2 1 — w^'^ 

- due to unstable filtrations with £ = 3, £^ = 1 and 1 < ^ 3: 

]:hi^o{z,Tf. 


Subtracting the terms above from (14. 2 p gives (15. lip . Subtracting further ^hi^Q{3z,3T) = 
3 6>i(62 3T)rj(3T) fro^^ (I5.1ip provides integer invariants in agreement with the definition (13. 2p . 

□ 

The recursive procedure explained above can be solved, such that hr,ci{z,T; J^^i) can be 
directly expressed in terms of the Hr'{z,T] f) with r' < r, without computing first the 
hr',ci{z,T; Js^i), and moreover giving more compact expressions. The solution follows from 
Ref. |32] (the solution to the recursion for vector bundles over Riemann surfaces) and Eq. 
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( 13. 6p one obtains: 



z — / ^ J. J. 

(ri,cii)A h{r£,C]^_^) = {r,cj), 1 = 1 

A'i'-^e,l>A'i+l'-^e,l 

/I \m-l ™ 

(5.15) = ^ ^ ^^;^2E.<, nr,(a,-a.) -Q^^^^^^^. J) 

(i-X.ajJH |-('-^,£i^) = (r,C2-C), i=l 



This becomes after carrying out the sums over [32 



-1) 



»Tl— 1 



m 

(ri ,ai)H Krm, am ) = {'", a) 

/ ^2M(ri,...,r,;a./n) 



where 



^-1 

(5.17) M(ri, . . . , r,; A) = ^(r, + r,+i) {(n + ■ ■ ■ + r,)A}, 



with {A} := A - [AJ. 

One can verify that Eq. fl5.16p for r = 3 is in agreement with Eqs. fl5.10p and (15. lip . As 
an example we give here h^fli^jT; Je,i)'- 



T 1 _|_ 7/;8 1 _|_ 7/;8 

/i4,o(^,r) = H,{T,z-f) + --—^H,{T,z-ff + -—-^H,{T,z-f)H,{T,z-f) 

z 1 — 1 — w'^ 



1 — w^^ 1 1 — w^^ 
(5.18) +- J- ^ //i(r, z; /)2 //^(r, /) + j-^ifi(r, z; /)^ 
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which is to be compared with: 

/i4,o(^,r) = H^{t,z; f) - [ — — + 



;i - (1 - U;12)2 2 (1 - (1 - «;24) 

11 11 f M 

'31^^ " 41^^ + 24 J 

2(1 + U'20) ^ 2 



(1 - ^16) (1 - U724) (1 - ^x;12) (1 - u;24) I _ yj 

+ 7:]hi,o{z,Tfh2fl{z,T;Js^i) 



24 



1-^16 2^ 

2(^10 + ^30) 2^18 ,,,,,,, 

+ 71 i:m71 ::5JT ^i.ol^,^") h2j{z,T;Je,i) 



(1 - (1 - u;24) (1 - m;12) _ ^24) 

— -^ + 1) hi,o{z,T)h3^o{z,T;J^^i) 

1 — W^^ J 

2 (w^ + v^^) 



\ — w 



24 



^1,0(2:, r) h3fl{z,T] Je,l). 



5.2. Wall-crossing. This subsection explains how to compute /i^.d {z^ t; J) for generic choice 
of polarization J from the generating functions for J = J^.i. The BPS invariants Q{T, w; J) 
for J differ in general from those for J = i, since sheaves might become semi-stable or 
unstable by changing the polarization. The change of the BPS invariants depends on the 
choice of Hirzebruch surface E^. Knowing how hr^ci{z,T; J) varies in C(Si) is particularly 
important for the computation of t; P^) since the blow-up formula is to be applied 

for the polarization Ji q = (f)*H, where H is the hyperplane class of (^see the next section). 
The change of the invariants can be obtained recursively from Eq. (15.81) after determining 
which filtrations change from semi-stable to unstable or vice versa. 
More quantitatively one has for J and J' sufficiently close: 

1 ^ 
AniT,w;J^J') = Yl lA ,(iri ..[ ^"^"^"""^^"^""'^'""'n^^r-'^"^) 

pj(ri)>rp,/(ri+i) 

(5.19) - |Aut({r.};.or ''^-""'"""''"n"'^'^'"-^''- 

Pj:(ri)hPj,(ri+i), 
pj(ri)^Pj(ri+i) 
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with |Aut({rj}; J) I defined below Eq. (15. 8p . Note that the invariants are evaluated in the 
chamber where the filtration is unstable. This makes this formula an recursive formula as it 
requires knowledge of Q{ri,w; J'), but since we are only interested in small rank this is not 
a serious obstacle. The recursion can be solved in various ways, e.g. [T2| [T3| [T6] or the Higgs 
branch analysis of Ref. p2] based on Ref. [25] . 

Since generating functions capturing wall-crossing are already described in the literature, 
the explicit expressions of hr^ai^, Jm,n) for r = 2 and 3, are presented here without further 
details. We have for r = 2 E]: 

i Ea,6ez 1 ( sgn((26 - (3)n - (2a - a)m) - sgn((26 -(3)- (2a - a)E) ) 

X (^y-{^-2)(2f)-/3)-2(2a-a) _ ^(£_2)(2b-/3)+2(2a-a)^ ^| (2(,-/3)2 + i(2fe-/3)(2a-Q) ^ ^^|2 ^ 

and for r = 3 [211120]: 

h3,l3C-af{z, T; Jm,n) = hs^/3C-af{z, T] J^^l) + 

\ ( sgn((36 - 2/3)r2 - (3a - 2a)m) - sgn((36 - 2/3) - (3a - 2a)e) ) 

^ j|y^-(i!-2)(36-2/3)-2(3a-2a) _ ^(£-2)(36-2/3)+2(3a-2Q) j ^^(36_2/3)2 + |(36-2/3)(3a-2a) 
xh2^bC'-af{z, r; S^, J|3fe-2/3|,|3a-2Q:|) hi^o{z,T). 

6. BPS INVARIANTS OF 

The Hirzebruch surface Si can be obtained as a blow-up : Si — ^ of the projective 
plane P^. Interestingly, we can compute the BPS invariants of P^ from those of Si from the 
blow-up formula. This formula is a remarkable result which states that ratio of generating 
functions of BPS invariants of a surface S and its blow-up S* is a (theta) function independent 
of 5 or J [301 El HE]- The underlying reason for this relation is that every semi-stable sheaf 
on S can be obtained from one on S by an elementary transformation along the exceptional 
divisor of the blow-up. 

Two subtle issues of the blow-up formula are (Proposition 3.4 of [30j): 

- the stability condition is /i-stability rather than Gieseker stability, 

- it involves the virtual Poincare functions X(r, w; J) of the moduli stack. 

To take these two issues into account let ^'^(r, w; J) be the invariant enumerating /i-semi- 
stable sheaves which is obtained from r2^(r, w; J) by addition of the Gieseker unstable sheaves 



BPS INVARIANTS OF SEMI-STABLE SHEAVES ON RATIONAL SURFACES 19 

which are /i-semi-stable using Eq. (15. 8p . Moreover, let X^(r,w; J) be the corresponding 
virtual Poincare function with corresponding generating function H^^^{z, r; S, J). 
The blow-up formula now reads [30j : 

Proposition 6.1. Let S be a smooth projective surface and (p : S S the blow-up at a non- 
singular point, with Ci the exceptional divisor of (p. The generating functions H^^_^{z, r; S, J) 
and r; S, J) are related by the "blow-up formula": 

Krc-kcM ^' ^*J) = BrA^^ r) H.^^Jz, T- S, J), 

with 



k 



The two relevant cases for this article are r = 2, 3: 

E„n^„ ,n „m'^+n'^+mn„ Am+2n 

[b.l) B2,k{Z, T) = -^-p , Bs,kiz, t) = -^-p . 

Note that Bj-A^i does not depend on S or J. 

The computation of hr^cii^, t; P^) from /ir,(^*ci-fcc(-2, t; Si) in general involves the following 
three steps: 

(1) Compute <^*ci-fcc;(-^' '^i "^i.o) by adding to hr^^*ci-kc{z,T\ Ji ^) terms due to sheaves 
on Si which are not Gieseker stable for Ji ,., but /i-semistable for (f)*H = Ji^q, and 
consequently compute ^-t^_^_j^A^i'^'^ Jifi) by adding the terms prescribed by Eq. 
(13. 5p . The generating functions and the factorial factors in Eq. (13. 5p combine these 
two steps very naturally into one. 

(2) Divide by 3,.^^, t) to obtain r; P^). 

(3) Determine /ir ci(-2^, P^) from r; P^) by reversing step (1). 

For Ci = pC+f, /3 = or 1, and J = Ji q, /i-stability is equivalent to Gieseker stability, and 
therefore steps 1) and 3) become trivial. For example, one can compute h^^niz, t; P^) either 
starting from h^j{z, r; Ji q) (which requires Conj ect ur e 14 . 1 1 and Eq. (15. 8p ) or h^^c+fi^^ 'T'-, Jifi)- 
One can verify that the first terms of both g-expansions of hs^niz, t; P^) are equal, which is in 
agreement with Proposition [6]T] A proof of the equality of these expressions for hs^niz, t; P^) 
would imply a proof of Conjecture 14.11 for (r, Ci) = (3,/) since h^j{z,T\ Js,i) is related to 
hsfi'+f{z, t; Ji^e) by the blow-up formula and wall-crossing. 
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When /i- and Gieseker stability are not equivalent, steps 1) and 3) are not trivial. We will 
first explain them for r = 2 following [29j. One obtains: 



Proposition 6.2. 



1 



b<0 
b=-l mod 2 



-/ii,o(^,r;p2)2. 



Proof. The only extended HN-filtrations which are Gieseker unstable for J = Ji ^ and fi- 
semi-stable for J = Ji o have i = C.^ = 2. For the parametrization Ci{E2) = hC — a/, the set 
of sheaves which is unstable for Ji ^ but yU-semistable for Ji q corresponds to 6 < and a = 0. 
This gives the second term inside the brackets. Consequently, step (2) divides by i?2,i(-2,T), 
and step (3) subtracts the yU-semi-stable sheaves which are not Gieseker semi-stable with 
£ = 2 and £^ = 1. □ 

Alternatively, one can compute h2fi{z,T]W'^) starting from h2fl{z,T; Ji^e). In that case the 
term due to step (1) in the brackets is ( b<o w^q^^^ + \] hi q{z, r)^, and one divides by 

V 6=0 mod 2 ^/ 

B2,o{z,t). Addition of ^hi^Q{z,T;¥^) provides the expected integer invariants, in agreement 
with [29]. Accidentily, the terms due to step (1) and step (3) can simply be incorporated by 
replacing Ji^^ by Ji,o in h2,i3c{z, r; Ji,e), and can be written in terms of the Lerch sum j2]. 

The remainder of this section discusses with r = 3. In terms of h^^ci^j 'T'-, Ji,e)i ^3,0(2^5 P^) 
is given by: 



Proposition 6.3. 



/i3,o(^,r;P' 



(6.2) 



1 



Bz,i{z,t) 
( 



( 



E 



w qi 



1 6<0 

y b=-2,—4 mod ( 



hi,o{z,T) h2fi{z,T; Ji^J 



+ 



E 



b 

1(7 g 12° 



1 6<0 

\ b= — 1 , — 5 mod 6 



+ 



E 



+ 



Ik 

w qi 



, fci,fe2<0 
\ fc2 — fc]^+l mod 3 



fc<0, 
A; — — 1,-2 mod 3 



hlfl{z,Tf 



-i/ii,o(^, r- P2)3 - ^ hi^^{z, T- P2) /i2,o(^, r; P^). 
o 2 
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The desired integer invariants are obtained from h^fii^z, r; P^) after subtraction of 0(82;, 3r 
5 e{6z 3t) • "^^^ ^^^^ non- vanishing coefficients are presented in Table [TlR They are in agree- 
ment with the expected dimension of A^(r) (12. 3p . 

C2 60 &2 &4 &6 ^8 ^10 bi2 bu &16 ^18 ho ^22 ^24 ^26 ^28 X 

3 1 1 2 2 2 2 18 

4 1 2 5 9 15 19 22 23 24 216 

5 1 2 6 12 25 43 70 98 125 142 154 156 1512 

6 1 2 6 13 28 53 99 165 264 383 515 631 723 774 795 8109 
Table 1. The Betti numbers bn (with n < dimcA^) and the Euler number 
X of the moduh spaces of semi-stable sheaves on with r = 3, Ci = 0, and 

3 < C2 < 6. 



Proof. The terms added to h^ ci^-, t; Ji^e) in the brackets are due to step (1). The last term 
on the first line and the term on the second line are due to filtrations with i = = 2. If 
one chooses Ci{E2) = bC — af as for r = 2, the set of sheaves which are unstable for Ji ,, but 
/i-semistable for Ji q corresponds to 6 < and a = 0. Similarly, the first term in parentheses 
on the third line is due to £ = = 3, and the second term due to i = 3 and £^ = 2. The 
sum of the terms in the bracket is Ht^^{z, r; Ji,e), and is divided by B^^ii^z, r) following step 
(2). Finally, step (3) corresponds to the last line. □ 

In a similar manner, /13 0(2;, r; P^) can be computed from h^fl[z,T\ Ji,^). Then the terms 
due to step (1) are for £ = 2: 



(6.3) 



and for 
(6.4) 



3: 



+ 2 A^'''Ki,o(^,r)/i2,o(2;,r;Ji,,) 



;j<o 

b— mod I 



+ 



E 



b<0 
fa— — 3 mod 6 



J 



E 



w 



1 fci,fe2<0 
\ k-^—k2 mod 3 



fc<0 
fc — mod 3 



^Note that these numbers are different from Refs. | [T71 [37], where Euler numbers of moduli spaces of 
sheaves and vector bundles are computed using localisation methods. This difference might be due to 
different choices of stability conditions. 
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